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ABSTRACT 

Fully sampled degree-scale maps of the ^^CO 2-1 and CO 4-3 transitions 
toward three members of the Lupus Molecular Cloud Complex — Lupus 1, 111, 
and IV — trace the column density and temperature of the molecular gas. Com- 
parison with IR extinction maps from the c2d project requires most of the gas 
to have a temperature of 8-10 K. Estimates of the cloud mass from ^^CO emis- 
sion are roughly consistent with most previous estimates, while the line widths 
are higher, around 2 kms^^. CO 4-3 emission is found throughout Lupus I, 
indicating widespread dense gas, and toward Lupus 111 and IV. Enhanced line 
widths at the NW end and along the edge of the B228 ridge in Lupus 1, and 
a coherent velocity gradient across the ridge, are consistent with interaction be- 
tween the molecular cloud and an expanding H I shell from the Upper-Scorpius 
subgroup of the Sco-Cen OB Association. Lupus III is dominated by the effects 
of two HAe/Be stars, and shows no sign of external influence. Slightly warmer 
gas around the core of Lupus IV and a low hne width suggest heating by the 
Upper-Centaurus-Lupus subgroup of Sco-Cen, without the effects of an H I shell. 
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Introduction 



The Lupus s tar forming; region, rece ntly reviewed by lComeronI ( 120081 ). lies about 150 pc 
from the Earth (ILombardi et al.l l2008al ). in the Gould Belt. It is immediately visible by 
inspection of optical photographs of the sk y, comprising a set of largely filamentary dark 
clouds. Based on optica l extinction maps (ICambresvl 119971 ). the 'Cores to Disks' Spitzer 
Legacy Programme (c2d, lEvans et al.lboosl . booj . I2OO9I ). has produced an d analyzed infrared 
images of three areas: Lupus I, III, and IV0. The c2d data products ( Evans et al. 2007 ) 
include: shorter-wavelength Infrare d Array Camera (IRAC) maps, used to identify and 



classify young stellar objects (YSOs. iMerm et al.ll2008l ): far-IR MIPS maps (IChapman et al. 



20071 ) tracing the thermal emission of the dust component of the mo 
extinction maps, based on the Two Micron All Sky Survey (2MASS; 



ecular clouds; and IR 



Chapman et al.l 120071 ) 



and on a combination of 2MASS and IRAC data (lEvans et al.l 120071 ). In this paper, we 
present maps of the c2d fields in Lupus, tracing the molecular hydrogen by the rotational 
emission of carbon monoxide (CO) and isotopically substituted ^^CO. 

The Lupu s comp lex has been extensively surveyed in the various 1-0 transitions of CO: 
Murphy et al.l (119861 ) first mapped the Lupus molecular clouds using the 1.2 m Columbia 



telescope with 0.5° resolution to show the extent of the molecular gas in Lupus and suggest a 
total mass of a few 10^ Mq. Improved maps in several transition s were obtained by NANTEN 
(2.7' resolution, but routinely undersampled) : iTachihara et al.l (1200 ll ) published ^^CO maps 
of the whole complex, mostly sampled at 8' spacing, with some areas at 4' spacing; ^'^CO 1-0 
maps (Tachihara et al.lll996l ) cover Lupus I and III with 8' spacing, while C^^O 1-0 maps 
( iHara et al.lll999l ) cover all the clouds in the complex at 2' spacing. iMoreira fc YunI (|2002l ) 
mapped about 200 arcmin^ of Lupus IV at sub-arcminute resolution and full-beam spacing 
in the 1-0 transitions of CO, ^^CO, and C^^O. The ^^CO 2-1 and CO 4-3 maps presented 
here are the first large-scale fully sampled low- J CO maps of Lupus, and the first mid- J CO 
maps of any kind. 

The large-scale str ucture of the Lupus c louds can be traced by near-IR extinction 
( ILombardi et al.ll2008bl ) and CO 1-0 emission (ITachihara et al.ll200ll ). The complex covers 
some 20° of galactic latitude (about 50 pc). At low latitudes {b < 10°), a large mass of diffuse 
gas contains denser filamentary clouds; at higher latitudes, the dense clouds (Lupus I and II) 
are more clearly separated. T he evolution of the Lup us clouds may have been driven by the 
influence of nearby OB stars (jTachihara et al.ll200ll ). The Lupus and Ophiuchus molecular 
clouds face one another across the Upper-Scorpius subgroup of the S cor plus- Cent aurus OB 
Association; the Upper-Centaurus-Lupus subgroup lies on the opposite side of the Lupus 



^Some authors use arable numerals: Lupus 1, 3 and 4 
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clouds from Upper-Sco. The H I shell around Upper-Sco, blown by stellar winds and a pre- 
sumed supernova 1.5 Myr ago fide Geuslll992l ). borders the NE side of the Lupus clouds; on 
the plane of the sky, the ridge which dominates Lupus I (B228, see below) lies just on the 
trailing edge of the H I shell. The much older U pper-Cen-Lup subgroup has driven an H I 
shell far beyond the Lupus clouds (Ide Geuslll992l ): it should have passed through the Lupus 
complex some 4-7 Myr ago, r oughly consistent with the ages of T Tauri stars in Lupus III 
and IV jMoreira fc Yunll2002h . 



Lupus I (Figure[I]) is dominated by B 228 (lBarnardlll927l ) or GF 19 (ISchneider &: Elmegreen 



19791 ). a long ridge running NW-SE, extending over about 2° (5 pc) parallel to the edge of 
the Upper-Sco H I shell. The molecular material appears to fall off steeply toward the center 
of the shell (at the NE edge), but there is extensive material on the other side of the ridge. 
Lupu s III (Figure [2]) is a long (about 4x1 pc) E-W cloud (GF 21. ISchneider fc Elmegreen 
19791 ) at the edge of the low-latitude cloud mass: at its eastern end, it curves up to the NE 
and breaks up into clumps; toward the west lie two embedded (but optically visible) Herbig 
Ae/Be stars, HR5999 (A 5-7) and the B6 HR6000 JComeronlboOsf ). iTachihara et~an J2002h 
consider Lupus III to be a cluster-forming structure with a dense head and spread-out tail, 
analogous to the p Oph and Cha I clouds. A couple of pc north of Lupus III, there is an- 
other smaller cloud, which appears to be a separate condensation in the low-latitude gas; 
we henceforth refer to this northern structure as Lupus IIIN. Lup us IV (Figure [3l) is the 



head of a filamentary struct ure running approximately E-W (GF 17. ISchneider fc Elmegreen 



19791 : iMoreira fc Yunll2002l ). comprising a small core, about 0.5 pc long (E-W), lying within 
a more diffuse extended cloud, ab out a pc across. The fila ment connects Lupus IV to the 
low-latitude gas mass to the East; ISchneider fc ElmegreenI describe it as a "chain of small 
faint globules". 

Individual dark clouds or cloud peaks have been identified by inspection of optical sur 



veys on scales ranging from degrees (e.g., B228) to arcminutes ( 


Sandavist & LindroosI 


1976|; 


Feitzineer & Stiiwe 


1984; 


Hartlev et al. 


1986; 


Bourke et al. 


1995a: 


Andreazza & Vilas-Boas 


1996; 


Lee & MversI 


I999I; 


Vilas-Boas et al. 


2000): these clouds are listed in Table [H Near-IR 



extinction maps of Lupus III with a resolution of 30" ( iTeixeira et al.ll2005l ) have been used 
to identify dark cores in more detail, which were then fitted with Bonnor-Ebert profiles. 
Single-pointing s pectra in several mole c ular lines have been observed toward extinction- 
selected clumps: ISandgyist fc LindroosI (19761) observed H2CO along four lines of sight at 
6 cm (6.6' beamwidth); iBourke et al.l (Il995br) searched fo r ammonia, but did not detect it, 
toward two globules (BHR120 & 140): Ivilas-Boas et al.l J200oh o bserved ^=^C0 and C^^O 
1-0, but their 0.8' b eamwidth is poorly- matched to ours; iLee et al.l (120041 ) observed a few of 
their Lupus I cores JLee fc MversI Il999h in CS 3-2 a nd DCO+ 2-1 (0 .7' beamwidth). Pre- 



viously identified CO cores, both those observed by IVilas-Boas et al.l . and those identified 
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from their C^^O maps by iHara et al.l (119991 ) . are listed in Table [2l The positions of these 
extinction cores and C^^O cores with respect to optical and ^^CO 2-1 emission are shown in 
Figures [TH51 Of the m embers of the c2 d small dark cloud sample mapped by AST/RO in 
^^CO 2-1 and CO 4-3 (ILohr et al.ll2007l ). three lie within the Lupus region, but outside our 
maps. 



Lupus I contains a few YSOs, concentrated to the B228 ridge (iMerm et al.l 120081 ). 
Lupus III contains a very dense cluster of YSOs near the He rbig Ae/Be stars, composed 



mainly of T Tauri st ars flAUen et al. 



object LupusSMMS (ITachihara et al 



200<i : iMerm et al.ll2008h. bu t including the Class 



20071 ). IComeron et al. I fl2009h identified a population 



of cool stars and brown dwarfs toward Lupus I and III, which are not strongly concentrated 
toward the mo lecular clouds. Th e Lupus IV region contains about as many YSOs as the 
Lupus I region (IMerm et al.ll2008l ). but they are older (Class II/III) and largely found outside 
the molecular cloud, with the exception of one flat-s pectrum YSO candid ate close to the 
cloud center. Even the dispersed population is absent (IComeron et al.ll2009l ). L upus IV thus 



appears to have less (or no) ongoing star formation, compared to Lupus I and III. lMerm et al. 



suggest that the Class II and III YSOs surrounding it represent an earlier generation of star 
formation (possibly associated with the passage of the Upper-Cen-Lup H I shell mentioned 
above), and that the very dense Lupus IV core is poised to form new stars. The y ounger 
YSOs (Class 0, I, and fiat-spectrum sources from the c2d samples, iMerm et al.ll2008l ) within 
the boundaries of our ^^CO 2-1 maps are listed in Table [H and plotted in Figures [THS 



Observations 



Lupus I, III, and IV were observed with the Antarctic Submillimeter Telescope and Re- 
mote Observatory (AST/RO. IStark et al.ll200ll ). during the period 2005 March-November, 
in the 461 GHz CO 4-3 and 220 GHz ^^CO 2-1 transitions. The telescope, receiver, and 



spectrometer systems are described by I St ark et al.l (120011 ): ^^CO 2-1 observations used the 
230 GHz SIS receiver, whose local oscillator system has been upgraded to use a synthe- 
sizer, frequency multiplied (xl9) by a Millitech multiplie r. CO 4-3 observations used the 
450-495 GHz SIS waveguide receiver (IWalker et al.lll992l ). Low- resolution acousto-optical 
spectrometers (AOSs, 1 GHz bandwidth, 0.7 MHz channel width) were used for all observa- 
tions. 

The large-scale maps presented in this paper were built up from multiple smaller maps 
('submaps'), mosaicked together as follows: each submap was centered on a point of a large- 
scale grid superimposed on the molecular cloud, the grid points separated by 24' and 12' 
for the 2-1 and 4-3 observations respectively. The submaps measure 28' x 29' {a x 6) and 
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12' X 12' respectively, giving substantial overlap between neighbouring 2-1 maps and some 
overlap (one row or column of spectra) between neighbouring 4-3 maps. 

The submaps use cell sizes of 1' and 0.5', significantly oversampling the telescope beam 
(3.3' at 220 GHz, 1.7' at 461 GHz). The spectra are position switched, taken in batches of 
4 or 5, all of which share a reference measurement, about a degree away. The 4 or 5 source 
spectra are all at the same declination (and hence at the same elevation, since AST/RO 
lies at the South Geographic Pole), but the reference position is the same for all spectra 
in a submap, and may thus differ in elevation from the source spectra by as much as a few 
arcminutes. The use of one reference position for each submap, however, allows the reference 
position to be checked for emission. Each row of each submap is built up from a number 
of contiguous blocks of reference-sharing spectra. Reference sharing increases the observing 
efficiency almost to that of on-the-fly mapping, but also shares one of its weaknesses: the 
sharing of a reference spectrum produces correlated noise in neighboring spectra, which can 
show up as artifacts in the map. Rather than the long stripes characteristic of on-the-fly 
mapping, correlated noise is manifested in the AST/RO maps as short (4- or 5-cell) horizontal 
blocks. 



3. Data Reduction 



The AST/RO observing system produces spectra calibrated onto the scale (IStark et al 



200ll ): because of the unobstructed off-axis optical design of the telescope, this is essentially 
equivalent to T^. Further reduction and processing, were carried out (using COMB, the stan- 
dard AST/RO data reduction progranj^, and PDljfl). The spectra have velocity resolutions of 
0.9kms~^ (^'^CO 2-1) and 0.4kms~^ (CO 4-3); they were checked for frequency shifts, and 
corrections were applied (corrections ranged from to almost 5kms~^ — see Appendix |A|) . 
Linear baselines were subtracted from all spectra, along with polynomial baselines if neces- 
sary: a few percent of the 4-3 spectra required non-linear baseline subtraction, and almost 
none of the 2-1 spectra (Table H]). 

Long-integration spectra were taken toward all reference positions, generally showing 
upper limits of ~ O.IK (^■^CO 2-1) and ~ 0.3 K (CO 4-3). Some 2-1 reference positions 
(for one submap in Lupus I and one in Lupus HI) contained substantial emission (~ IK), 
and these reference spectra were added back into all the spectra in the relevant submap. 



^http://www. astro. umd.edu/^mpound/comb 
■^the Perl Data Language, http://pdl.perl.org 
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Because of the overlap between submaps, there are over a thousand pointings on the 
sky toward which two or more spectra were measured at different times. These duphcate 
pointings can be used to estimate the internal consistency of the flux calibration at the 
telescope. Distributions of the ratios of peak between spectra with the same pointings 
are shown in Figure H] for both 4-3 and 2-1 data. The typical inconsistencies between these 
spectra (estimated by the standard deviation of the distributions of ratios) are about 
30% for CO 4-3 data, and about 20% for ^^CO 2-1. 

3.1. Map Making 

Initial map making was carried out by nearest-neighbor sampling the submap data onto 
one large grid, with overlapping observations co-added. Each pixel of these maps corresponds 
directly to the spectrum observed by the telescope; without smoothing, inconsistencies be- 
tween neighbouring submaps show up as clear edges, which are much more obvious in the 
4-3 than in the 2-1 data. 

The overlap between submaps allowed these inconsistencies in the 4-3 data to be checked 
and corrected. Individual submaps with significantly different flux scales to their surround- 
ings were identified and rescaled. These rescalings were checked visually: a submap was only 
rescaled if it was obviously higher or lower than its neighbours, if the measured scaling was 
consistent with the visual one, and if the rescaling improved the appearance of the map when 
it was applied. Rescaling corrections of < 10% did not significantly improve the appearance 
of the map, and were not implemented. Fewer than half the submaps needed correction, 
mostly by factors of 0.7-1.4 (in line with the estimate of internal consistency given above), 
and the largest correction factor was ~ 1.8. No corrections needed to be applied to the 
2-1 data. The need to rescale probably arises from imperfect estimation of the atmospheric 
transparency in the calibration procedure, which works better when the atmosphere is fairly 
optically thin (e.g. at 220 GHz) than at higher opacities (e.g., at 461 GHz). 

The corrected 4-3 data (and uncorrected 2-1 data) were then gridded (as above) to 
produce final unsmoothed data cube^, which were used to generate the maps in this paper. 
All the quantities mapped are calculated directly from the data, without profile-fitting. For 
example, the line width is estimated by the ratio of integrated intensity to peak (which 
overestimates the FWHM by about 10% if applied to a perfect gaussian line). Sample spectra 
from the data cubes are shown in Figures [5H71 



publicly available in FITS format at http://www.astro.ex.ac.uk/people/nfht/resourcesJupus.html 
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Integrated intensity maps (Figures [81 4T01) . channel maps (integrated over 2 kms ^ chan- 
nels, Fig ures [TT]-[T3|). and peak maps (Figures [T^1 - [T6|) were all produced by bessel smooth- 



mg (e.g. iJenness Sz Lightfootll2000l ). which gives the 'most fair' representation of the sky as 
observed by a single-dish telescope. However, if significant emission extends to the edge of 
the map, this technique tends to produce a roll-off of emission in the outermost pixels of the 
map, and this effect is visible in some of the Lupus maps. The effective telescope diameter 
was picked to balance resolution against noise: 1.5 m for the 2-1 maps, and 1.2 m for the 
4-3, yielding overall beam sizes of 3.8' and 2.3', respectively. Maps of the velocity centroid 
(Figures [T71 - [T9|) and line width (Figures [20] - [22l) use both unsmoothed data (in grayscale) 
and Gaussian-smoothed data (contours), with Gaussian FWHMs of 5' (2-1) and 5.75' (4-3) 
yielding overall beams of 6' in both transitions. 

The nominal pointing accuracy of AST/RO is about 1', but comparison of the ^^CO 2-1 
map of Lupus IV with the optical extinction structure (see Figure [3]) suggests a systematic 
offset in de clination. Discre pancies of similar magnitude can be seen in other AST/RO data 
from 2005 (ILohr et al.ll2007l ). The Lupus IV maps have therefore all been shifted southwards 
by 1.5'. The submillimeter map is well correlated with the millimeter- wave one, and so it 
has been shifted as well. There is no evidence of similar systematic shifts in the Lupus I or 
Lupus III data. The publicly available data cubes do not have this shift applied to them. 



4. Results 



4.1. Overview of the Lupus Clouds 



Vilas-Boas et al.l (120001 ) estimated C^^O 1-0 optical depths in the range 0.1-0.5 toward 
extinction peaks, implying that ^^CO 2-1 will no longer be optically thin i n dense parts of 



the Lu pus clouds. Maps of CO 2-1 peak follow the extinction maps of iChapman et al. 



(120071 ) better than maps of integrated intensity do, suggesting widespread breakdown of 
the optically thin approximation. The maximum peak toward Lupus I and IV is about 
3.5 K in each map, implying a minimum excitation temperature, T^, of 8K; the peak in 
Lupus III reaches 5 K near the bright stars, so > 10 K in this area. iTachihara et al.l (119961 ) 
adopted 13 K for their analysis of ^^CO 1-0 toward Lupus I. A comparison of the lower ^^CO 
contours to the most recent c2d extinction maps suggests that excitation temperatures of 8- 
10 K will reproduce the extinction-traced column density quite well. While higher excitation 
temperatures are not ruled out by the molecular hne data, they would lead to lower column 
density estimates which would be less compatible with estimat es from extinction, wh ich is 
probably the least biased estimator of the true column density (iGoodman et al.l 120091 ). 
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Gas masses are estimated from the ^^CO 2-1 data as follows: the ratio of the peak 
to the maximum possible brightness temperature (for assumed T^) yields an estimate of 
the peak optical depth Tq^^-*; this is multiplied by the line width to estimate / rdv, which 
is proportional to the column density of ^^CO {N^^^\ Appendix [B|) . A^(^3) is converted to 
a total gas column density by assuming conversions between A^(^^) and extinction Ay (see 
below), and betAv een Ay and the column density of hydrogen (iVn = 1.37 x lO^^Aycm"^, 
Evans et aliboooh . Assuming a distance of 150 pc, the total gas mass enclosed in a 1 arcmin^ 
pixel is then calculated. 

The relationship between ^'^CO column density and extinction can be expressed as 
jY(i3) ^ ^(13) y. xo^^(y4y — Avfi) cm~^, with the extinction threshold Ay^ reflecting the 
minimum column density required for the presence of ^fCO. The most commonly used pa- 



Frerking et al. 



cm 



mag 



rameters come from a study of p O ph bvl. , 

Avfl = 1.6 mag; other studies (e.g. Bachiller fc Cernicharol 1986 : Lada et al. 1994 ) have found 



values of x^^^^ in the range 2.2-2.7, and ex tinction thresholds of 0-1.6 mag. Recently, a com 



prehensive study of the Perseus complex ( Pineda et al.ll2008l ) found x*^^^^ = 2.4 cm ^ mag 



Av,o = 1.7mag overall, and x*-^^-* = 1.9 cm~^ mag~^, Ayfl = 1.2 mag toward the West End of 
the complex (PWE), comprising several dark clouds, which is likely to be the best analog to 
the Lupus clouds. 

Estimated masses of the Lupus clouds above Ay thresholds of 2 and 3 mag are given in 
Table along with the mean and dispersion of the line widths of material with Ay > 2. (The 
line width distribution for Ay > 3 is very similar, varying only by about 0.1 kms~^.) Masses 
have been calculated for = 8, 9, 10 K and for three ^^CO-to- Ay relations: PWE, Perseus, 
and FLW. The quoted mass is the estimate using = 9K and PWE, and the positive 
and negative errors give the full range of mass estimates. Other estimates of mass and line 
width from the literature are also tabulated. Previous isotopically substitute d CO results 
relied, directly or indirectly, on the older ratio of Nh to Ay (1.87 x 10^^ cm~^, iBohlin et al. 



19781 ). and have been rescaled to use the same ratio as in this work. The extinction mass of 
Lu pus III has also been r escal ed to assume a dis tance of 150 pc instead of the 200 pc used 
by 



Chapman et all (120071 ) and Evans et al.l (l2009f ). 



The masses of Lupus I, III (including Lupus IIIN) and IV are estimated to be 280, 
150, and 80 Mq respectively (for Ay > 2). The extinction mass of Lupus I (251 Mq) is 
similar to the ^^CO 2-1 estimate; those of Lupus III and IV are significantly greater (250 and 
120 Mq), but still consistent with the ^^CO 2-1 error range. This suggests that the Lupus 



comp lex has spatial variations in the ^^CO-to- Ay ratio, as Perseus does ( Goodman et al. 
2009 h. By contrast, the ^^CO 1-0 mass estimate toward Lupus I is 880 Mq (iTachihara et al.l 
1996,) . The ^'^CO 1-0 map covers a much larger area than ours (almost 10 deg^) and includes 
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significant emission to the SE and SW of our map boundaries. Witfiin tlie confines of tlie 
^^CO 2-1 map, we estimate tliat tlie ^^CO 1-0 mass sliould be 600-700 Mq, still about 
50 % more than th e uppe r limit of our mass range. A similar discrepancy was reported 
by iMoreira fc Yun whose mass estimate for Lupus 11 is about a factor of 3 lower 

than that of iTachihara et al.l : but their mass estimate for Lupus IV is quite similar to ours. 
Possible reasons for the discrepancies include the undersampling of the NANTEN maps 
(which will tend to ov erestimate the are a of structures) and differences in the calibration 
of the two telescopes. iHara et al.l (119991 ) measured only 240 toward Lupus 1 in C^'^O 
1-0, which is much closer to our estimate. The additional mass in the ^^CO 1-0 map might 
comprise tenuous gas in which the 2-1 transition is not fully excited, and whose column 
density is low enough for it not to be detectable in C^'^O 1-0. The ^^CO 1-0 mass estimate 
toward Lupus 111 (220 Mq) is consistent with our data and close to the extinction result, 
while the C^^O 1-0 estimates toward Lupus 111 and lllN (80 and 5 Mq) are much lower. 

The average line widths for the Lupus clouds (Tableland Evans et al. 2009)f| have been 
deconvolved by subtracting the AOS channel width in quadrature, and, at 1.5-2.7 kms~\ 
are broader than the ^^CO 1-0 widths repor ted at emiss i on pe aks by ITachihara et al.l (119961 ). 
much broader than the 0.9 kms~^ given by iHara et al.l (119991 ) as the typical line width, and 
also broader than the C^^O 1-0 line widths toward individual cores (Table [2]). C^^O 1-0 is 
expected to display narrower lines than ^^CO 2-1, since the more optically thin transition 
arises mainly from the denser gas in the center of the cores, which generally has a lower 
velocity dispersion. It is harder to see why there should be a discrepancy between ^^CO 1-0 
and 2-1 line widths. However, the 1-0 measurements are taken from two spectra extracted 
from the map, so the data are insufficient to investigate the disagreement. Despite these 
differences, the various line width measurements follow the same pattern: Lupus 1 and 111 
have quite similar line widths. Lupus 111 N has broader lines and Lupus IV narrower. 

Based on the line widths of ^^CO 2-1 and the velocity gradients observed toward the 
clouds, it is clear that they are gravitationally unbound (although their general filamentary 
shape makes the standard formula for virial mass hard to apply). Lupus IV, however, may 
be marginally bound due to its high column density, narrow lines and low velocity gradients: 
the lower limit to its virial mass is estimated to be ~140 Mf7) , while the upper limit to the 
mass derived from the ^^CO 2-1 map is ~15OM0. iHara et al.l (|l999l ) found almost all of the 
dense cores that they identified in C^^O 1-0 to be gravitationally unbound; the greater line 
widths in ^'^CO 2-1 reinforce this conclusion. 



^The tabulated line widths are slightly different from those in the earlier paper: 1 Evans et al.l averaged 
the mean and median line widths and took the difference between mean and median to be the error range; 
we only report the mean. 
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4.1.1. CO 4-3 Emission 

Because of their small areal coverage, concentrating on ^'^CO 2-1 emission peaks, the 
CO 4-3 maps do not provide the same extensive overview of the clouds. The 4-3 map of 
Lupus I, however, covers most of the B 228 ridge. Throughout most of the ridge, the peak 
lies between 1 and 2 K, corresponding to excitation temperatures of 7-9 K (assuming that 
the CO emission is optically thick throughout, with a filling factor of unity). The ends of the 
ridge may be warmer, with peak of 3-4 K implying excitation temperatures of 10-12 K. 

The 4-3 maps of Lupus III and III N show the majority of the gas to be similar to the 
ridge in Lupus I, with ~ 8 — 10 K. Within the main Lupus III cloud there are two maxima: 
one with peak of 4K [T^ ~ 12 K), while the other peaks at 9K. This latter one appears 
very compact, and so its excitation temperature may be significantly higher than the 18 K 
minimum implied by its peak T^. The 4-3 map of Lupus IV covers only the very center of 
the cloud, with peak of 3-4 K (T, ~ 10-12 K). 



4.2. Lupus I 



Towards the B 228 ridge, maps of extinction (I Chapman et al. I l2007h . peak ( Figure [1] 
and integrated intensity (Figure [8]) ar e not well correlated. The extinction map is the most 
robust estimator of column density (iGoodman et al.ll2009l ): it shows a very strong peak 
in the middle of the ridge (DC 338.8+16.5), slightly less to the SE (339.2+16.1), and a 
broad low plateau to the NW (338.7+17.5). The ^^CO 2-1 peak map has similar peaks 
toward all three dark clouds, while the CO 4-3 peak map shows no significant peak 
toward 338.8+16.5 in the center. The integrated intensity map is very different from the 
extinction map, with the strongest emission in the NW followed by 339.2+16.1 in the SE, 
and 338.8+16.5 comparatively weak in the center (hardly appearing in the CO 4-3 integrated 
intensity map). 

Previous molecula r line studies haye found a velocity gradient o f about 0.3 kms^^pc"^ 



along the B228 ridge (ITachihara et al.l l200ll : IVilas-Boas et al.l 120001 ). but the channel and 
velocity maps (Figures [TT] and [T7|) show very little evidence of a smooth velocity gradient 
along the ridge. Rather, the two large clumps at the NW and SE end of the ridge have 
velocities that differ by about 3 kms~^, while the ridge between them lies at an intermediate 
velocity. On the other hand, a velocity gradient across the center of the ridge shows up 
very clearly in both 2-1 and 4-3 maps: the gradient of about 1 kms~^ pc~^ is quite strong, 
compared not only to previously reported gradie nts in this cloud, but als o to the gradients 



of 0.1-0.4 kms pc found in areas of Taurus (iMurphy &: MyersI Il985l ). The gradient is 
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smooth and coherent over the whole of this region — an extent of more than 1 pc — and 
appears quite hnear, which would be cons istent with sohd-body rotation around the long 
axis of the filament (iGoodman et al.lll993l ). A coherent pattern in this region may also be 
seen in the line width maps (Figure [20]) . which show an increase in line width toward the 
leading (NE) edge of the ridge. There is no sign of a velocity gradient across the NW half 
of B228, which suggests that the NW part of the ridge may be distinct from the center. 

One of the peaks in the ce nter of the ridge, D C 3388+165-5, is associated with IRAS 
15398-3359: this Class YSO ( Shirley et al. 200oh has a cornpact j^V) molecular outflow 
with a dynamical timescale of about 2kyr ( Tachihara et al.l Il996l ). It lies about 0.5 pc 
behind the steep edge of the ridge, in the E-W elongated core that forms the southern half 
of DC 338.8+16.5, and coincides with local maxima in integrated CO 4-3 emission (Figure[H]), 
in peak of both transitions (Fig fT^ . and in CO 4-3 line width (Figure [211]) . 

At the SE end of the ridge, the NE spur of DC 339.2+16.1 is blueshifted from the 
ridge, while the rest of the clump is at about the sam e velocity, su g gestin g that the spur 
might be a separate cloud superposed on the ridge. iLee fc Myerd (119991 ) identified four 
optical extinction peaks in this clum p: DC 3392+16 1-1 to -4. They associated Peak 1 with 
IRAS 15420-3408/HT Lup, a CTTS JComer6n|[2008l ). which is optically visible as a nebulous 
patch between p eaks 1, 2 and 3 in FiglH The other 3 'starless' peaks were observed in CS 3-2 
and DCO+ 2-1 Lee et al.l[20oi ): all three have T^(CS)~ 0.5 K; DCO+ was only detected 
(with T^(DCO^)~ 0.2 K) toward peaks 3 and 4. DCO"*" is a high-density tracer, and its 
presence toward peaks 3 and 4 along the ridgeline, and absence from peak 2 at the start of 
the NE extension, suggest that, even in this complex structure, dense gas is concentrated 
to the ridge. The CO 4-3 data, with better spatial and velocity resolution, show additional 
structure in this clump: a patch of strong emission close to DC 3392+161-4, about 10' across, 
is blueshifted by about 1 kms~^ with respect to the main ridge structure, leaving a cavity in 
the emission at ambient velocity. Molecular outflow is a possible explanation, but since there 
is no obvious red lobe, this would require the Y SO to lie b e hind the cloud. Blue asymmetries 
in molecular lines can be caused by infall, but iLee et al.l (120041 ) searched unsuccessfully for 
signs of infall toward DC 3392+161-4. 

The parsec-scale ring-shaped structure to the SW of B228 is not covered by the CO 
4-3 map. However, its ^^CO 2-1 emission looks quite similar to that seen in the bulk of the 
ridge, and the detection of DCO"*" toward the dense core embedded in the SE edge of the ring 
suggests a significant amount of dense gas. The structure may also have a strong velocity 
gradient: the channel maps (Figure [TT]) show two complementary semicircular structures in 
adjacent channels, i.e. separated by about 2kms~\ equivalent to a velocity gradient of at 
least 1 kms"^ pc~^. In the line width map (Figure [20]) . the redder NW side of the structure 
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has somewhat broader hnes. iLee et aLl (12004 ) observed an extinction core at the SE edge of 
the ring: they found no signs of infall, but both CS 3-2 and DCO"*" 2- 1 are stronger in this 
core than in any of the cores observed in B228. iTeixeira et al.l ( 120051 ) found a ring-shaped 
structure in extinction in Lupus III, with a diameter of about 5', which they interpret as the 
remnants of the molecular cloud that formed the nearby cluster of young stars. By contrast, 
this structure is larger (20') and is not associated with a known cluster. 



4.3. Lupus III 



CO emission toward Lupus III in both transitions is dominated by a co mpact source 
in the west of the cloud, close to core F identified by ITeixeira et al.l (120051 ). Proceeding 
from E to W are: Core F; the compact CO source (slightly blueshifted from the bulk of the 
cloud); the HAe/Be stars HR5999 and 6000; and further ^^CO 2-1 emission extending E-W. 
This latter ernission shows no particular sign of the ring structure found in extinction by 
Teixeira et all J2005h . probably because the central void is smaller than the 220 GHz beam. 



There is a significant velocity gradient across this area, of about 2 kms^^ pc^^, getting redder 
to the west, which is not simply due to the blueshifted emission of the compact source to 
the east. 



To the east of the compact source, extinction cores A-E (ITeixeira et al.l |2005|) blend 
together at this resolution: the integrated intensity map (Fig [9]) shows two E-W spurs off 
the emission peak, the northern containing cores C and D, and the southern containing E. 
Further east still, the northern spur continues to core B, while A lies just NE of DC 3397+92- 
3. The strongest H13qo+ emission in this area (~ 1 Kkms ^) is found toward Core E; 
Core D, with about h alf the H^^CO"^ emission, contains the Class YSO Lupus 3 MMS 
dTachihara et al.lboOTl l 



Lupus IIIN, the small (< 0.5 pc) cloud core lying a few pc N of Lupus III, shows 
surprisingly complex structure: ^'^CO 2-1 channel and velocity maps (Figures [T^] and [THl) 
show the core to be lying in a larger (~ 2 x 1 pc) structure, elongated SE-NW, with a 
significant velocity change over the long axis. The velocity gradient is dominated by a 
sudden change just SE of the cloud core, also seen as a sharp edge in the reddest 2-1 channel 
map. 
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4.4. Lupus IV 



The E-W elongated core of Lupus IV has 3 peaks along its length: the central and 
western ones correspond to the extinction peak which reaches Ay ~ 24 mag. The peaks are 
most obvious in the ^^CO 2-1 channel maps (Figure [T3|) : the peak map (Figure [T6|) shows 
them with much lower contrast than the extincti on maps as the ^ ^ CO t ransition is optically 
thick {Av ~ 24 mag corresponds to r^^^^ ~ 8-10). iMoreira fc YunI J2002h found the ^^CO 1-0 
line to be optically thick as well, saturating with respect to 100 fim IR emission. The eastern 
peak has broader lines yielding similar integrated i ntensity fFigure [TO|l: at finer velocity 
resolution, these lines are found to be double peaked ( iMoreira &: Yurul2002l ). Although there 
is quite strong CO 4-3 emission throughout the core, and the velocity structure of the 4-3 
emission is similar to that of ^^CO 2-1, the integrated 4-3 emission is concentrated toward 
the E and W condensations only. Two elongated structures running NE of the core are also 
visible in the CO 4-3 peak map. 

The strong velocity gradients seen in Lupus I and III are absent from Lupus IV. Some 
velocity structure is evident: the central condensation in the core and the two elongated 
NE structures a re blu eshifted (by 1-2 kms^^) compared to the E and W condensations. 
Moreira fc YunI (120021 ) reported significant velocity structure in the E-W direction, but our 
velocity maps show that this E-W structure is not found throughout the cloud. 



4.5. Comparison of CO transitions 
4.5.1. Whole clouds 

CO 4-3 emission from the Lupus clouds can be compared to that of ^^CO 2-1 by 
Gaussian smoothing the CO 4-3 map and plotting the integrated intensities (Figure [23]) and 
peak (Figure against one another, point-by-point. The peak comparison is noisier 
than that of integrated intensity, and is only shown for Lupus I and III, which have stronger 
emission. 

The large 4-3 map of Lupus I provides an overview of the bulk of the molecular gas. Clus- 
ters of points close to the origins of the integrated intensity and peak plots reflect the noise 
levels in the maps. The majority of the data are somewhat correlated, occupying regions of 
the plots bounded by (4-3)/(2-l)> 0.5 (both plots), / T^(4 - 3)dv < 6 K km s'^ (integrated 
intensity), and T^(4 — 3) < 2.5 K (peak). The 4-3 emission is effectively saturated, while 
the ^'^CO 2-1 emission takes a wide range of values. This suggests that even regions of quite 
low column density contain sufficient dense gas to emit optically-thick CO 4-3 radiation (the 
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critical density of CO 4-3 is of order 10^ cm ^, but see Sec. 15. II) . 

Another population of molecular gas can be identified toward Lupus I: its ^^CO 2-1 
peak is similar to that of the bulk population, while its CO 4-3 peak is greater 
(T^(4 — 3) < 4K). This suggests that the gas is warmer, but not denser than the bulk 
of the cloud. The integrated intensities of the transitions, however, are well correlated; 
this probably arises from increased line widths in both transitions, and thus from increased 
velocity dispersion in the gas. Inspection of the maps of integrated intensity, peak and 
line width toward Lupus I (Figures [HI [IH [201) shows that this component, with enhanced 
integrated intensity, peak CO 4-3 and line width, is found in the NW end of B228. 

The much smaller CO 4-3 maps toward Lupus III, IIIN, and IV provide fewer data. 
However, two components can still be identified in Lupus III: one component consistent 
with the Lupus I bulk population, and the other with stronger emission. In contrast to its 
counterpart in Lupus I, both integrated intensity and peak of the brighter component are 
well correlated. This brighter component constitutes the compact peak close to the HAe/Be 
stars in Lupus III. Heating by the nearby stars could increase the peak and hence the 
integrated intensities of both transitions, without a significant change in the line width. The 
molecular gas in Lupus HI N and IV occupies a similar region of the integrated intensity plot 
to the bulk population of Lupus I, suggesting broadly similar physical conditions. 

Some positions show significant CO 4-3 emission without ^'^CO 2-1: the spectra toward 
one of these positions are shown in Figure [3 and the integrated intensity and peak plots 
show that this occurs in all the clouds. The combination of strong CO 4-3 emission and 
weak or absent ^^CO 2-1 implies warm gas with low column density; since the CO 4-3 map 
areas were chosen to cover ^^CO 2-1 peaks, there may be more of this population outside 
our 4-3 maps. 



4-5.2. Dense Cores 



Spectra toward specific cores have been nieasur ed by NANTEN and the Swedish-ESO 
Submillimeter Telescope (SEST): [Hara et al.l (119991 ) list the C^'^0 1-0 lin e parameters for 



the sp ectrum at the peak of each dense core identified from their map, and [Vilas-Boas et al. 



( 12000[ ) list parameters for both ^^CO and C^^O transitions toward extinction-selected cores. 
Because the AST/RO maps are fully sampled, it is possible to obtain equivalent spectra 
(albeit with a beamwidth of 3.3', compared to the 2.8' NANTEN beam, and SEST's 0.8' 
beam) and compare the peak (Figure [25l) and integrated intensities (Figure [261) . 



The peak in CO 4-3 and ^"^CO 2-1 toward most of the cores lie in the same region 
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of the plot as the bulk component of molecular gas described above. Four cores are clearly 
brighter than the rest: the two identified in Lupus III are both associated with the bright 
compact source, the one in Lupus I lies at the SE end of B 228, and the one in Lupus IV is in 
the middle of the central core. The majority of the extinction cores have ^^CO 2-1/1-0 peak 
ratios in the range 0.3-0.7, inconsistent with the 0.7-2 expected from LTE (Appendix [B|) . 
This could be caused by beam-dilution compared to the sub-arcminute SEST beam (but see 
below). The exception (one of the Lupus III cores associated with the bright compact source) 
has a ratio of almost 3: it may be large enough not to suffer beam dilution, and warm 
enough to be close to the high-temperature line ratio limit. While the Lupus I and III cores 
are well mixed in the plot, the Lupus IV cores have rather low line ratios. There is no 
particular reason for them to be more beam diluted, so the lower line ratio may reflect a 
lower temperature. The plot of ^^CO 2-1 peak against that of C^^O 1-0 does not show 
any clear separation between the extinction and C^^O core samples. This is surprising, since 
the C^^O cores observed by NANTEN should suffer almost as much beam dilution as the 
AST/RO data, compared to the extinction-selected cores observed by SEST. The line ratios 
of the cores lie mainly between unity and ~ 3; correcting for the factor-of-2 discrepancy 
estimated above, this sug gests line ratios of abou t 2-6, in line with LTE estimates for C^^O 
optical depths of 0.1-0.5 Jvilas-Boas et al.lboooh . 

The CO 4-3 and ^'^CO 2-1 integrated intensities of the cores are more scattered, but 
the high-T^ cores in Lupus III also have the highest integrated intensities. In contrast, the 
high-T^ cores in Lupus I and IV have more average integrated intensity, suggesting that the 
increased (and hence temperature and/or density) is not accompanied by increased line 
width. The two more cores in Lupus I with large ^^CO 2-1 integrated intensity also lie in 
the SE of B228. The ratios of the two ^^CO transitions (2-1/1-0, for the extinction-selected 
cores only) lie between about 0.7 and 1.5 (with a few between 1.5 and 3), which is consistent 
with our expectations from LTE, but not with the peak ratios above. Line parameter s 
were estimated from the SEST data by fitting Gaussian line profiles ( IVilas-Boas et al.ll2000l ). 
in contrast to the AST/RO and NANTEN data, for which integrated intensity and peak 
were measured directly, and the line width estimated by the ratio. This difference in analysis 
could cause a systematic discrepancy, with the SEST analysis yielding lower and broader 
line widths. The cores with higher ratios (1.5-3) are all in Lupus I or III, but are scattered 
throughout the clouds. Ratios of ^^CO 2-1 to C^^O 1-0 integrated intensity are around 
3-6 for the C^^O-selected cores observed with NANTEN (consistent with the peak ratios 
above), but the same ratios toward extinction cores are generally higher. These high ratios 
generally arise from rather low C^^O 1-0 integrated intensities (< 0.5 Kkms~^). 
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4.6. CO Emission toward YSOs 



The YSO populat ion of Lupus is dominated by Class II and III objects (I Comer on 



20081 : iMerin et al.l 120081). I put younger objects are also found there. The Class 0/I/F sources 
identified by iMerin et al.l lying within the ^^CO 2-1 maps are listed in Table [3], together 
with the 2-1 line parameters toward them and, where a pplicable, CO 4-3 pa rameters; two 
objects that have been identified as background galaxies (IComeron et al.ll2009l ) are excluded. 
Spectra toward the YSOs are shown in Figure [271 all 2-1 spectra are pointed within 0.6' 
of the YSO position, and all 4-3 spectra are within 0.25'. Spectra toward some YSOs do 
not show any significant emission: values of below ~ IK (in either transition) should be 
treated as noise. 



Of the 4 YSOs in Lupus III with very low or non-existent 2-1 emission, IComeron et al. 



( 120091 ) estimate 3 to have ages of order 100 Myr. The rest of the YSOs have line widths 
similar to those of the surrounding molecular gas. This does not rule out the presence of 
line wings due to outflow: the Class source Lupus 3 MMS has quite average line widths in 
both transitions, but the spectra themselves clearly show wings. The fiat-spectrum source 
J154506. 3-341738 (in Lupus I) h as broader -than-average lines; nebulosity prevented it from 
being measured in the optical by IComeron . 



5. Discussion 



5.1. The Bulk Molecular Material 



The estimates of column density, and hence mass, derived above assume d a blanket 
temperature of 9 ± 1 K rather than the 13 K ass umed bvlTachihara et al.l (119961 ) for Lupus I 
and the 12-17 K estimated toward Lupus IV by Moreira fc Yun ( 20021 ). both based on opti- 
cally thick CO 1-0 emission. 8 K is the minimum excitation temperature consistent with the 
peak ^^CO 2-1 seen toward the majority of the gas in Lupus I, III, and IV, and with the 
peak CO 4-3 over most of the B 228 ridge. Adopting excitation temperatures close to the 
minimum implies assumptions of optical thickness, LTE, and a filling factor close to unity, 
but is required for consistency with the extinction maps (lEvans et al.ll2007l ): the = 2.5 
and T^(^^CO 2-1)= 1.5 K contours are quite similar. This implies T^; < lOK for both PWE 
and FLW ^^CO-to-Ay relations, although it does allow a higher toward Lupus I for the 
Perseus ^^CO-to-Ay parameters. 



The gas temperature estimate is consistent with estimates for many dark clouds. IVilas-Boas et al 



(I2OOOI ). using ^^CO and C^^O 1-0 transitions, estimated the excitation temperatures of their 
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sample of dense cores in Lupu s to lie in the range 7 -15 K, with the majority of good esti- 
mates being colder than 10 K; I Clemens et al.l (Il99ll ) estimated gas temperatures toward a 
large sample of small dark clouds from CO 2-1, and found that the majority wer e colder than 



10 K. Models of cold dark clouds yield dust temperatures <10K at the center (lEvans et al. 



200ll : IZucconi et al.ll200ll ). While the dust temperature is higher toward the cloud surface, 
gas temperatures are l ower than dust temperatures for low column and volume densities (e.g.. 
Doty fc Neufeldl 119971 ) ■ so the ^^CO 2-1 and CO 4-3 transitions need not be dominated by 
the cloud centers. The CO 1-0 estimates mentioned above are likely to be dominated by the 
cloud surfaces, which may be warmer than those of the majority of dark clouds due to ex- 
ternal heating by the nearby OB association. The peak CO 4-3 toward Lupus IV implies 
a Tx of 10-12 K, compared to ^ 8K from ^^CO 2-1 . This may reflect a combination of 
external heating, as suggested by lMoreira &: YunI (120021 ). with high enough density to couple 
the gas and dust temperatures more effectively than in Lupus I or III. 

The column density of H2 through out most of the B 228 ridge is within a factor of 2 
of 5 X 10^^ cm~^ (jChapman et al.l 120071 ). The width of the ridge on the sky is about 10', 
or 0.4 pc; if B 228 is assumed to be a filamentary cloud, with a similar depth, the average 
volume density n is a few 10^ cm~^. This is close to the critical density ric of the 2-1 
transitions of CO and its isotopologues, supporting the assumption of LTE used throughout 
this work, but an order of ma g nitud e lower than the critical density of CO 4-3 (a few 



10"^ cm ^, Appendix [B]) . lEvand (119991 ) found that significant emission in many transitions 



could arise from gas with volume density more than an order of magnitude lower than Uc 
(although CO was not included in that study), but the similarity between the excitation 
temperatures of ^^CO 2-1 and CO 4-3 suggests thermalized emission, and thus a volume 
density close to Uc- The line emission is pervasive, which argues against its arising from cores 
or a high-density center of the ridge. However, if the bulk of the ridge material were taken 
to be close to the critical density of CO 4-3, the depth of the ridge along the line of sight 
would be an order of magnitude lower than its width in the plane of the sky, which seems 
implausible. It is more likely that the CO 4-3 emission arises from a small fraction of dense 
gas found in clumps throughout the molecular cloud; or from a thin warm shell around the 
outside of the cloud; or that subthermal emission from the bulk of the molecular gas can 
account for the CO 4-3 lines seen throughout the cloud. The latter possibility, in particular, 
may be checked by modelling the emission. 



5.2. Cloud Cores 



Vilas-Boas et al.l (120001 ) estimated the optical depth of C O 1-0 toward their dark cloud 
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sample to lie in the range 0.1-0.5 by comparing the ^^CO and C^^O lines, and assuming an 
abundance ratio of 5.5. A large fraction of their cores had line ratios in excess of 5.5, 
inconsistent wit h their assumed i sotopologue ratio; however, ratios up to 8 are possible 
(Appendix [B]). iHara et al.l (Il999l ) estimated somewhat lower optical depths toward their 
C^^O cores, assuming = 13 K; at T-^ = 9 K , their estimated optical depth range becomes 
quite consistent with that of IVilas-Boas et al.l . The ^^CO 1-0 optical depth toward the cloud 
cores is likely to range from 0.5 up to 2-4, depending on the abundance ratio. 



Teixeira et al.l (120051 ) estimated volume densities of a few 10"^ cm ^ for the dense cores 



they identified in Lupus III. If such densities are widespread among the dense cores in Lupus, 
the CO 4-3 transition is close to LTE, and it becomes possible to estimate the expected value 
of (T^) 4,-3 /{T^) 2-1 toward the cores. At = 9K, the high-column-density limit of the ratio 
(both transitions optically thick) is ~0.4; as the optical depth of ^^CO becomes moderate, 
the ratio increases toward unity. However, at low optical depths, CO 4-3 may no longer 
be thermalized, making this estimate invalid. Most of the cloud cores have 4-3/2-1 ratios 
clustered around 0.5 (Figure [25]), with a few more at unity or above, and one over 1.5. It 
is difficult to achieve a line ratio significantly above unity under LTE: even at low optical 
depth, the temperature would need to be at least 20 K. 



5.3. Lupus I 

The central condensation in the B228 ridge (338.8-1-16.5) has high extinction, moderate 
peak and integrated in ^^CO 2-1, and CO 4-3 emission quite similar to the bulk of 
the cloud. The line width is also quite similar to that of the bulk cloud, so the enhancement 
in ^^CO 2-1 is largely due to the increased peak T^. Elevated gas temperature would likely 
show up in the CO 4-3 emission, so the increased is largely due to the increased column 
density as the ^^CO 2-1 transition becomes optically thick. This dense cloud seems to have 
more in common with the bulk material around it than it does with the emission peak to the 
NW, having similar temperature and line width. It is also associated with two embedded 
YSOs, including a known outflow source. 

The strong, coherent velocity gradient seen north and east of 338.8-1-16.5 runs perpen- 
dicular to the leading edge of the ridge (i.e., in the direction of the H I shell's expansion); 
the leading edge (i.e., toward the center of the shell) has a greater line width than the gas 
behind it. It is difficult to rule out temperature and optical-depth effects combining to mimic 
a velocity gradient, but similar patterns are seen in both ^^CO 2-1 and CO 4-3 maps: a truly 
optically thin tracer is required to confirm the gradient. If the apparent gradient truly refiects 
the kinematics of the gas, the total change in velocity across the ridge is about lkms~^, or 
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half the hne width. However, the change in velocity across a Jeans length (around 0.1 pc for 
these clouds) is only about O.lkms"^, which is unlikely to add significant support against 
collapse. Indeed, 338.8+16.5 seems to be part of this velocity field, and contains a known 
Class YSO. 

The integrated intensity maximum at the NW end of B 228 is due to the combination of 
enhanced peak and broader lines. Throughout this area, the peak CO 4-3 is > 3 K, 
implying a gas temperature of at least 10 K. Although the 4-3 emission probably does not 
sample the whole of t he gas column, there i s supporting evidence for a warmer temperature: 
the Spitzer IR maps (IChapman et al.l 120071 ) show th is end of the filarn ent to be bluer than 
the central part, with strong 24 /^m emission (see also lMerfn et al.ll2008l ). which could be due 
to a higher dust temperature . No Class 0/I/F YS Os are found in the northern part of the 
ridge, and only one Class III (IChapman et al.ll2007l ). so there are no obvious internal heating 
sources. The nearby Upper-Sco OB subgroup lies to the NE of Lupus I, but is unlikely to 
be causing the heating, since no similar effect is to be found in the center of the ridge. 

At the SE end, the picture is complex: the extinction and peak for both transitions 
are high at the SE end of 339.2+16.1, while line width and integrated intensity peak further 
NW. The strong enhancement of peak CO 4-3 at the SE end suggests increased tem- 
perature as well as column density. The associated YSOs all lie to the NW in this area, so 
internal heating seems unlikely. 



5.3.1. Interaction with the Upper-Sco Shell 



The H I emission from the Upper-Sco shell lies at a similar velocity (3-9 km s ^. Ide Geus 



19921 ) to the Lupus clouds, consistent with their being in contact with one another, and 
Lupus I being dynamically affected by the shell. If such an interaction is going on, the 
combination of continuity and conservation of momentum require that the sum of pressure 
and momentum flux is conserved by the interaction. The limitations of both molecular and 
atomic data make it impossible to look for clear diagnostics of interaction between the shell 
and ridge, but some rough estimates can be made. For both phases. 



(P + pv^) /k 



n 



cm '^K Vcm \ Vkm s ^/ Vkm s ^ 



+ 



where = 170 for H I, and 340 for H2. Based on the H I data jde Geuslll992l ). this quantity 
may be estimated to be of order 10^ cm~^K, about 90% of wh ich is the mom entum flux 
component due to the expansion velocity of the shell (10kms~^, Ide Geuslll992l ). P + pv for 
the molecular gas in Lupus I is also around 10® cm^'^K, but is approximately evenly divided 
between pressure and momentum flux (due to velocity gradients). 
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The rough similarity in this sum between the H I and molecular gas is consistent with 
interaction between them. If the H I shell is indeed affecting Lupus I, it is likely doing so 
through the momentum flux of its expansion. This transfer of momentum could be causing 
the velocity gradient across the B228 ridge (in the direction of the H I she ll's expansion) ; 
this process could be analogous to the 'streamers' in the Ophiuchus complex (Ide Geuslll992l ). 
The enhanced line width at the NW end of B 228 and along the leading edge of the ridge 
(implying increased pressure in the molecular gas) is also consistent with the effect of the 
H I shell. More detailed studies of both H I and molecular gas may support the idea of 
Lupus I being affected by the Upper-Sco shell: better estimates of the volume density of H2 
and higher-resolution maps of H I are required. 



5.4. Lupus III 



The brightest part of Lupus III near HR5999 and 6000 was mapped by lTachihara et al. 

(120071 ) in the millimeter- wave continuum, C^^O 1-0 and H^^CO^ 1-0, all with higher reso 



lution than the AST/RO data. The C^'^O emission peak coincides with the CO peak in our 
maps, and extends to the east, with another E-W elongation a few arcminutes to th e north. 
The continuum map, together with a near-IR extinction map (INakajima et al.ll2000l ). shows 
the same southern E-W structure extending further east, while the northern structure breaks 
up into two clumps, the denser western one containing a Class protostar (Lupus 3 MMS). 
In the H^^CO"'" map, there is no emission at the CO peak; the emission peaks strongly to 
the east in the southern E-W structure, while the clumps in the northern structure show up 
as smaller peaks. 



Tachihara et al.l (120071 ) suggest that the lack of H^^CO^ emission at the CO peak can 



be accounted for by a long path length through gas with volume density significantly lower 
than the critical density (~ 10^ cm~^ for H^^CO"*"). They estimate a column density of 
1.2 X l O^^cm"^, equivalent t o Ay of about 13 mag, which is consistent with the extinction 



maps (IChapman et al.l 120071 ). The peak in CO 4-3 (about 9K) implies > 18 K, 
while the 6K peak of ^'^CO 2-1 only requires Tx > 11 K. Compared to the rest of the 
Lupus clouds, where the CO 4-3 and ^^CO 2-1 excitation temperatures are quite similar, 
this is a significant discrepancy. It could be explained by the CO peak being compact, as 
suggested by the CO 4-3 peak map, so that the ^'^CO 2-1 measurement is beam-diluted. 
Alternatively, the CO 4-3 could be tracing an outer shell heated by the HAe/Be stars. 
However, the ^^CO 2-1 emission is also optically thick (for = 11 K, Ay of 13 mag implies 
an optical depth of about 8), and so will trace similar material. These temperature estimates 
tend to support the argument that the gas is too warm for depletion of H^'^CO^ to explain 
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the lack of emission toward the CO peak (ITachihara et aLll2007l ). although the estimates are 
unlikely to apply to the center of the clump. The strong CO 4-3 emission implies that the 
transition is thermalized, so a significant amount of th e gas in the clump must have volume 
density ^ 3 x lO'^cm"^. While ITachihara et al.l (120071 ) derived an average volume density of 
lO^cm"'^, there must be enough significantly denser gas to thermalize the CO 4-3 line, but 
not enough to excite the H^'^CO''' 1-0 transition. 



In projection, Lupus III lies far away from any part of Sco-Cen (ITachihara et al.ll200ll ). 
in contrast to Lupus I and Lupus IV (discussed below), so the HAe/Be stars probably 
influence it far more than the OB associations. However, there is evid ence for Lupus III 



being further away than Lupus I and IV, possibly as much as 50 pc (e.g. lComer6nll2008l ). If 
this is the case, then Lupus III could lie behind Sco-Cen, and could be affected by either or 
both of Upper-Sco and Upper-Cen-Lup. The large line widths seen toward Lupus IIIN could 
be caused by external influence in the same way as the broader lines seen toward Lupus I. 



5.5. Lupus IV 



The peak in CO 4-3 toward Lupus IV of about 4K, implying > 12 K, occurs 
around the outside of the extinction peaks, which have peak around 3K (T^. > 10 K). 
This suggests that the outside of the clump is significantly externally heated. The ^^CO 
2-1 peak implies > 8 K, but even this measure is unlikely to sample the extinction 
maxima properly, since it will be dominated by the outer layers of the structure. The 
CO 4-3 temperature estimates are in line with those seen at the NW end of B 228, where 
Lupus I seems to be strongly affected by Upper-Sco. L upus IV is on the oppo site side of the 
Lupus complex to Lupus I, Upper-Sco and its H I shell (ITachihara et al.ll2001l). bu t faces the 



Moreira fc Yun 



(l2002f ) suggested 



Upper-Cen-Lup subgroup, which lies to the W and SW. 
that Lupus IV was shaped by the influences of both subgroups, and noted that some of 
the velocity gradients they saw in Lupus IV were along the vector toward Upper-Cen-Lup. 
Much of the enhancement in the peak CO 4-3 lies on the S and W sides of the extinction 
peaks, which is consistent with a picture of external heating by the radiation field from 
Upper-Cen-Lup. There is, howeve r, no nearby H I shell , the Upper-Cen-Lup shell having 
passed the Lupus clouds long ago (IMoreira fc Yurul2002l ). If interaction with the H I shell 
causes the enhanced line width at the NW end of B 228, the lack of any such interaction in 
Lupus IV would be consistent with its rather low line widths. 
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Conclusions 



Fully sampled degree-scale maps of the ^'^CO 2-1 emission toward the Lupus I, III, and 
IV clouds trace the column density and temperature of the gas, the transition becoming 
optically thick iii the c l oud cores. The pea k is well correlated with the near-IR extinction 



( IChapman et al.l 120071 : lEvans et al.l 120071 ). and a comparison of the two suggests that the 



Tachihara et al. 


1996; 


Moreira & Yun 


2002; 


Teixeira et al. 



20051 ). This estimate is fairly robust to changes in the relationship between ^^CO column 
density and Ay. Estimates of the cloud masses from the ^^CO maps are reasonably consistent 
with those derived from extinction mapping. The differences between these estimates vary 
greatly from cloud to cloud, and suggest that there may be si gnificant spatial variat ion in 



the ^^CO-to-y4y relationship, as found in the Perseus complex (iGoodman et al. 120091). The 
line widths of ^^CO 2-1 toward the clouds are higher than previous estimates (lEvans et al. 
20091 ). around 2kms~^, with Lupus IIIN rather broad and Lupus IV rather narrow. 



Fully sampled CO 4-3 maps covering most of Lupus I and small regions of Lupus III 
and IV trace dense gas: the peak generally indicates excitation temperatures quite close 
to those of ^^CO 2-1, and hence that the transition is largely thermalized. This suggests 
that the volume density > lO^cm"^, although modelling will be required to ascertain the 
required density. CO 4-3 emission is pervasive toward Lupus I (the map of which covers a 
large area), implying that this dense gas is found either throughout or all around the outside 
of the cloud, although it may comprise a fairly small fraction of the cloud mass. 

The physical conditions of the molecular gas vary along the B 228 ridge in Lupus I. At 
the NW end, the gas has broader lines and probably higher temperature than in the bulk 
of the cloud; the column density is not particularly high and there is only one Class III 
YSO. In t he center of the ridge , the dark cloud 33 8.8+16.5 is associated with recent star 
formation (^Tachihara et al.lll996l : IShirley et al. I l2000h : in this coherent velocity gradient 

of about 1 kms~^ pc~^ runs across the ridge. The SE end of the ridge is complex, with YSOs, 
enhanced line width and integrated intensity on the NW side, and column density (and 
possibly temperature) peaking to the SE. The enhanced line widths and velocity gradient in 
B 228 are consistent with a dynamica l interaction between Lupus I and th e H I shell around 
the Upper-Sco subgroup of Sco-Cen (Ide Geuslll992l : iTachihara et al.l 1200 ll ). 



To the north of Lupus III, the small cloud Lupus IIIN has similar characteristics to the 
bulk of the other clouds, albeit with broader lines and significant velocity structure. Lupus III 
itself contains a compact CO peak which is probably heated by the nearby HAe/Be stars 
HR5999 and 6000. The gas in this clump contains sufficiently dense gas to thermalize the 
CO 4-3 transition (n^ ~ 3 x lO'' cm~^), but not to thermalize the H^^CO^ transition mapped 
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by iTachihara et al.l (120071 ). which would require n > lO^cm"^. The rest of Lupus III seems 
to have quite similar physical conditions to those in the rest of the clouds, and shows no 
particular sign of being affected by the nearby OB subgroups. 



Lupus IV contains peaks of very high column density (jChapman et al.l 120071 ) associated 
with slightly warmer gas temperatures (10-12 K). These temperatures are estimated from the 
optically-thick CO 4-3 transition, which is strongest around the extinction cores, suggesting 
significant external heating. The average line width of ^'^CO 2-1, however, is significantly 
low er than those of Lupus I and III. Lupus IV faces the Upper-Cen-Lup subgroup of Sco-Cen 



and iMoreira &: Yud (120021 ) suggested that it is influenced by the OB stars; the Upper-Cen- 
Lup H I shell passed by the Lupus clouds a few Myr ago, so Lupus IV is more likely to be 
affected by the radiation field from the OB stars. 

Despite the basic similarities in their physical conditions, the three clouds have signifi- 
cant differences: Lupus I appears to be strongly affected by external thermal and dynamical 
influences from the nearby Upper-Sco OB association, and does not display widespread star 
formation. Lupus III shows no sign of external influence — parts of the cloud are heated in- 
ternally by its own young stars. Lupus III N seems entirely quiescent, yet has a large average 
line width. Lupus IV has the greatest column density and the narrowest average line width, 
has almost no star formation as yet, and may be heated externally by the Upper-Cen-Lup 
OB association. 

A detailed spatial comparison of CO and extinction maps will yield more accurate 
estimates of the physical conditions of the Lupus clouds, as well as mapping the variation in 
the ^^CO-to-Ay ratio. Mapping of additional CO transitions is crucial to the understanding 
of the clouds: more optically-thin lines (C^^O and even C^^O) are particularly important. 
Maps of the Upper-Sco H I shell with comparable resolution to the Lupus maps are required 
to look for more definitive signs of interaction between the shell and the molecular clouds. 
More sophisticated radiative transfer calculations are beyond the scope of this work, but are 
required to use the CO 4-3 emission as a proper constraint on the physical conditions of the 
gas, and hence to understand the structure of the Lupus clouds and how they are affected 
by the local environment. 
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A. Frequency Calibration 

Because the AOS backends are analog devices, laser mode hopping causes shifts in their 
frequency calibration. A number of such shifts occurred over the nine-month period during 
which these data were taken. 

The fundamental frequency calibration for each spectrometer was obtained once, by 
connecting a frequency synthesizer to the IF input of the AOS to obtain the AOS channel 
width and the channel number of a fiducial frequency. The frequency scale thus defined was 
used for all data in this paper, but some corrections had to be made. The change in channel 
width caused by a mode-hop is negligible, but the entire spectrum is offset by a few channels. 

The frequency shifts were tracked with a number of fiducials. Each AOS includes a 
comb generator which is usually used to obtain several frequency calibration scans per hour, 
giving excellent frequency tracking. However, the comb generator failed during 2005, so 

other frequency standards had to be used. For the CO 4-3 observations, the mesosphcric 
CO 4-3 line is so strong that it can be picked up without a switched measurement. So 
the 'sky' spectra, used to estimate the sky temperature for calibration, show the line at 
an antenna velocity close to zero. The mesosphcric CO abundance has a strong seasonal 
variation, and became so weak at the end of the austral winter (around September) that 
it could no longer be seen in the sky spectra. Finally, repeated spectra were taken toward 
the compact H II region NGC 3576. This source has velocity structure on the scale of the 
AST/RO beam, so pointing uncertainties translate into velocity uncertainties. These three 
fiducials were combined to track the frequency scale. The majority of ^^CO 2-1 spectra were 
corrected by l-3kms~^, and some were corrected by up to 5kms~^. About 90% of the CO 
4-3 spectra were corrected by < lkms~^, with the remainder corrected by 3.9kms~^. 

Frequency shifts also show up as inconsistencies in the maps. The 4-3 observations 
showed no obvious inconsistencies, but some shifts had to be applied to the ^'^CO 2-1 data: 
channel maps of Lupus I showed that data taken in 2005 November had a significant un- 
corrected frequency shift compared to earlier, better-calibrated data. In addition, the long- 
integration spectrum toward one of the reference positions with significant 2-1 emission was 
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shifted by about half a channel with respect to the map spectra into which it was added. 
This was not corrected, because the facility for combining spectra in COMB only handles 
integer-channel shifts, and the effect of the shift is negligible, even in the channel maps. 



B. CO Emission 

The CO and ^^CO data in this work are analyzed under the assumption of LTE (i.e., 
the excitation temperature, T^, and gas kinetic tempera ture, Tk, are the same). For the two 



transitions considered in this work. Equation (14.46) of iRohlfs fc WilsonI (119961 ) yields 

22.1 



1- (l + ri^) 



for CO 4-3 (neglecting the cosmic microwave background term), and 

^ 10.6 



lnfl+ 



r^(2-l)+0.21 



for ^^CO 2-1. These equations are correct if the transitions are optically thick, and com- 
pletely fill the beam; otherwise they underestimate T^. In Lupus, the gas is cold enough 
{Tx ^ hv/k) that is significantly different from T^. 

The column density in the lower level of the ^^CO 2-1 transition is given by 

n[^^^ = 9.69 X 10^^ — - / rdv cm-^ 

1 1 - exp(-10.6/T,) J 

and this can be converted to the total column density of ^'^CO by correcting for the partition 
function: N^^^'^ /n[^^^ ranges from 2.1 to 2.3 for of 7-10 K. 

In LTE, the ratio of ^^CO 2-1 to 1-0 emission depends only on temperature. If both 
transitions are optically thick, the ratio is simply the ratio of brightness temperatures at 
different frequencies for a given T^, and will range from 0.7 (at 7K) to unity (at high 
temperature). In the optically thin case, this ratio is multiplied by the ratio of optical 
depths, which also depends on via the partition function: this ratio ranges from 1 to 2 
between 6K and 15 K, with a high-temperature limit of 3. The 2-1/1-0 line ratio should 
therefore range from 0.7 to about 2 in Lupus. 

The ratio of ^^CO 2-1 optical depth to C^^O 1-0 optical depth is just the ratio of ^^CO 
2-1/1-0 optical depths, multiplie d by the ^^CO/C^^O ab undance ratio: the ^^C/^^C isotope 



ratio is 62 ± 4 in the local ISM flLanger fc Penziaslll993l ). but the double ratio ^^CO/C^^^O 
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is not as well known. Combining the solar isq/i ^q j-atio of 500 (Zinnerl Il996l) with the 
local ISM i^C/i^C ratio yields i^CO/C^^O ~ 8, but llanger fc PenziasI Jl99o[ ) point out that 
i2q^i3q g^j^^ i6q^i8q g]2ould track one another, being similarly d ependent on star formation 
history, and so a solar ratio (^^CO/C^^O ~ 5.5. lMyers et al.lll983h may be more appropriate. 
Thus the ratio of optical depths should fall in the range of 4 to about 16. The ratio of ^^CO 
2-1 emission to C^^O 1-0, however, is complicated by the fact that the C^^O transition is 
likely to be fairly optically thin, while the ^^CO transition will have moderate to high optical 
depth. This will tend to reduce the ratio: conditions in Lupus are likely to yield emission 
ratios as low as 2-4. 

Below a critical volume density rig, LTE fails (Tj, < Tk)- The critical density itself 
depends on physical conditions: the effective spontaneous emission rate Aij is reduced at 
high optical depth, yielding a lower ric, and collisional transition rates are temperature- 
dependent. In the optic ally thin limi t, nc(C04-3) varies from 2.9 x lO'' cm~^ at 40 K to 
4.2 X 10^ cm~^ at 10 K (jjanserull995l ): the cr itical densities of CO 1-0 and 2-1 transitions 
are a few hundred and a few thousand cm~^ ( iRohlfs fc Wilson! 1 19961 ). 
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Table 1. Dark clouds in Lupus 
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Table 1 — Continued 
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Table 2. Molecular Clouds in Lupus 
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1.2 


0.5 




Lu2 


15 39 56 


-34 42 50 








2.0 


3.0 


1.5 








4.6 


3.3 


0.7 


1.9 


1.0 


0.5 




Lu3 


15 40 10 


-33 40 07 


1.8 


4.5 


2.5 


1.8 


3.0 


1.6 








5.5 


9.9 


1.7 


0.5 


0.8 


1.6 


338.8+17.2 




15 40 14 


-33 38 41 


2.0 


5.0 


2.5 


2.4 


5.7 


2.4 


0.5 


0.8 


1.3 
















Lu4 


15 40 32 


-34 39 40 








2.3 


5.0 


2.2 








3.3 


2.9 


0.8 


1.2 


0.6 


0.5 




Lu5 


15 42 03 


-33 46 33 


1.9 


4.3 


2.2 


1.5 


4.0 


2.7 








3.0 


4.7 


1.5 


0.7 


0.5 


0.7 




Lu6 


15 42 04 


-33 50 36 


1.4 


2.5 


1.8 


2.1 


5.2 


2.4 








1.9 


1.9 


0.9 


0.5 


0.4 


0.8 




Lu7 


15 42 24 


-34 09 02 


1.8 


3.3 


1.8 


2.9 


5.3 


1.9 








4.0 


3.6 


0.9 


1.3 


1.1 


0.8 




Lu8 


15 42 35 


-33 52 50 


1.7 


4.6 


2.7 


2.5 


6.0 


2.4 








4.5 


6.2 


1.3 


1.1 


0.8 


0.7 


338.8+16.5 




15 42 35 


-34 08 58 


1.6 


2.8 


1.8 


2.8 


7.2 


2.6 


1.9 


2.2 


1.2 














339.0+16.7 




15 42 48 


-33 53 56 


1.5 


4.0 


2.7 


2.3 


6.6 


2.8 


1.4 


1.6 


1.1 
















B228 


15 43 02 


-34 09 06 


1.8 


3.3 


1.9 


3.0 


5.2 


1.7 








5.2 


5.9 


1.1 


0.8 


0.8 


1.0 




Lu9 


15 43 10 


-34 13 50 


1.6 


2.3 


1.5 


2.4 


4.4 


1.8 








4.2 


3.6 


0.8 


0.9 


0.4 


0.4 


339.1+16.1 




15 44 59 


-34 18 08 


2.2 


5.8 


2.6 


3.0 


9.1 


3.1 


1.0 


1.5 


1.1 
















LulO 


15 45 06 


-34 17 39 


2.2 


6.3 


2.8 


2.9 


8.6 


3.0 








3.4 


3.3 


0.9 


0.8 


0.5 


0.6 


339.1+15.9 




15 45 30 


-34 25 51 


3.2 


4.9 


1.5 


2.9 


5.9 


2.0 


2.3 


1.5 


0.7 














Lupus III 




lu30 


16 07 51 


-39 11 12 








2.4 


4.9 


2.1 








3.2 


2.2 


0.7 


0.3 


0.2 


0.7 


339.6+9.3 




16 08 53 


-39 06 26 


9.0 


15.5 


1.7 


5.6 


13.1 


2.3 


1.9 


2.4 


1.2 
















lu31 


16 09 08 


-39 03 55 


4.2 


8.0 


1.9 


3.3 


8.0 


2.4 








1.5 


1.8 


1.1 


0.4 


0.2 


0.4 




lu32 


16 10 19 


-39 12 16 


2.3 


3.0 


1.3 


2.2 


5.2 


2.3 








4.3 


3.9 


0.9 


0.3 


0.2 


0.7 




lu33 


16 10 27 


-39 05 18 


2.5 


3.7 


1.5 


2.1 


2.9 


1.4 








5.1 


4.1 


0.8 


1.2 


0.7 


0.6 




lu34 


16 11 16 


-39 02 38 








0.9 


1.9 


2.1 








1.0 


1.0 


1.0 


0.5 


0.2 


0.4 




lu36 


16 11 28 


-39 00 21 








1.3 


3.2 


2.5 








4.6 


3.7 


0.8 


1.2 


0.7 


0.5 




lu35 


16 11 36 


-39 04 18 








1.4 


2.7 


2.0 








4.4 


2.8 


0.6 


0.6 


0.5 


0.7 



Table 2 — Continued 



Cloud 




02000 


52000 




CO 4-3 




13 


CO 2- 


-1 




C^^O 1-0 


a 




i^CO 1- 


-Ob 




C^^O 1-0 


b 


Notes 


Hara 


VMF 








/ 


AV 




/ 


AV 


Tr 


/ 


Ay 




/ 


AV 




/ 


Ay 




340.7+9.7 




16 11 58 


-38 04 23 


1.5 


6.0 


4.1 


2.3 


5.4 


2.4 


0.5 


0.9 


1.8 














Lu III N 


Lupus IV 




lu24 


16 00 18 


-42 03 47 








1.4 


2.6 


1.8 








4.1 


2.4 


0.6 


0.8 


0.3 


0.4 






lu23 


16 00 57 


-42 04 55 








1.6 


2.9 


1.8 








4.2 


3.8 


0.9 


1.2 


0.7 


0.6 




336.4+8.2 




16 00 57 


-42 03 16 








2.1 


3.7 


1.7 


2.1 


1.4 


0.5 
















336.7+8.2 




16 01 46 


-41 52 34 


3.6 


5.4 


1.5 


3.4 


6.5 


1.9 


2.0 


2.0 


0.8 
















336.9+8.2 




16 02 34 


-41 41 54 


2.2 


2.9 


1.3 


2.0 


3.7 


1.8 


1.5 


1.3 


0.8 


















lu25 


16 02 36 


-41 42 26 


1.2 


1.2 


1.0 


1.9 


3.3 


1.7 








4.9 


2.5 


0.5 


0.6 


0.2 


0.4 




336.7+7.8 




16 03 12 


-42 07 43 








1.8 


3.7 


2.0 


0.7 


0.6 


0.9 
















336.8+7.9 




16 03 37 


-42 00 55 








1.5 


3.1 


2.1 


0.7 


0.7 


0.8 


















lu26 


16 04 10 


-42 00 40 








1.4 


2.7 


1.8 








5.2 


3.8 


0.7 


1.1 


0.6 


0.5 





HHara et al.] l ll999l) 



IVilas-Boas et al.l l|2000h 



00 
00 
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Table 3. Class 0/I/F YSOs in Lupus 



Designation 


Class'' 


O/NIRb aaooo 


<52000 




CO 4-3 




13 


CO 2- 


-1 


SSTc2d'' Other 










/ 


Ay 


T* 


/ 


AV 




Lupus I 


J154214.6-341026 


I 


Low 15 42 14.6 


-34 10 26 


1.8 


2.2 


1.2 


2.5 


4.7 


1.9 


J154301. 3-340915 IRAS 15398-3359 


0<^ 


15 43 01.3 


-34 09 15 


4.3 


8.2 


1.9 


2.8 


4.8 


1.7 


J154506. 3-341738 


F 


Neb 15 45 06.3 


-34 17 38 


2.1 


5.9 


2.8 


2.9 


8.6 


3.0 


Lupus III 


J160703. 9-391112 


F 


IR, Old 16 07 03.9 


-39 11 12 








0.4 


0.9 


2.4 


J160708.6-391407 


F 


V, Old 16 07 08.6 


-39 14 08 








0.6 


1.8 


3.2 


J160829.7-390311 Szl02 


I 


Note 16 08 29.7 


-39 03 11 


2.6 


3.4 


1.3 


2.7 


6.1 


2.2 


J160831. 1-385600 


F 


16 08 31.1 


-38 56 00 


0.6 


1.1 


1.7 


0.3 


0.2 


0.8 


J160846. 8-390207 


I 


16 08 46.8 


-39 02 07 


4.2 


7.8 


1.9 


2.7 


6.0 


2.3 


J160918. 1-390453 Lupus3MMS 


O'^ 


Neb 16 09 18.1 


-39 04 53 


3.5 


6.2 


1.8 


3.3 


8.0 


2.4 


J160934.2-391513 


F 


16 09 34.2 


-39 15 13 


2.3 


4.1 


1.8 


1.7 


3.8 


2.2 


J161013. 1-384617 


F 


IR, Old 16 10 13.1 


-38 46 17 








0.4 


0.6 


1.4 


J161027.4-390230 


F 


Low 16 10 27.4 


-39 02 30 


2.6 


4.4 


1.7 


1.9 


4.5 


2.4 


J161204.5-380959 


F 


V, Low 16 12 04.5 


-38 09 59 








1.2 


4.3 


3.7 


Lupus IV 


J160115. 6-415235 


F 


Neb 16 01 15.6 


-41 52 35 


4.5 


6.4 


1.4 


3.6 


6.2 


1.7 



HMen'n et al] ||2008| ') 

From lComeron et al ] l|2009l '): low, below main sequence; Neb, hidden by nebulosity; IR, suspected IR excess; old, age of order 
100 Myr; V, suspected veiling; note, discussed individually. 

'^see the text for details 
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Table 4. Areas mapped and quality of spectra for both transitions 



Region Area(deg^) Area/cells Spectra Non-linear^ % non- linear^ 



2-1 



Lupus I 


2.57 


9258 


11856 








Lupus III 


L87 


6740 


8412 


12 


0.1 


Lupus IV 


0.92 


3300 


3300 








Total 












Lupus 


5.36 


19298 


23568 


12 


0.05 






CO 4 


3 






Lupus I 


0.76 


10964 


12067 


939 


8 


Lupus III N 


0.04 


625 


650 


10 


2 


Lupus III 


0.17 


2401 


3240 


150 


5 


Lupus IV 


0.09 


1225 


1425 


15 


1 


Total 












Lupus 


1.06 


15215 


17382 


1114 


6 



^spectra requiring non-linear baseline subtraction 
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Table 5. Masses and Line Widths in Lupus 



Quantity 


Lu I 


Lu IIIN 


Lu III 


Lu III + 
Lu IIIN 


Lu IV 


M/Mq{Av > 2) 
M/Mq{Av > 3) 
AV (kms-i) 
rms AV^ (kms^^) 


1501^° 
2.3 
0.7 


2.7 
0.9 


2.2 
0.9 


1501^^° 

2.4 
0.9 


8011° 
301^° 
1.5 
0.4 


Previous Work 


M/M0 (i^CO 1-0)^ 880^ 
AV^/kms-H^^CO 1-0)^ L9 

m/Mq{^^co i-oy 

M/Mq (C^^O 1-0)'^ 240 
A^/kms-HC^^O l-O)'^ 0.6-1.7 
M/Mq{Av > 2Y 251 


5 

1.8 


1.2 


220 
1.7 

80 

250^ 


<150 

160 
0.5-0.9 

120 



Tachihara et al.l (119961 ) 



''Covering a larger area than the present work; see Section 14.11 



Moreira fc YunI fl2002f ) 



iHara et al.l (jl999|) 



Evans et al.l (120091 ) and N. J. Evans II (2009, private communication) 



^Assuming a distance of 150 pc 
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Right Ascension (JSOOD Oj 



Fig. 1. — Optical (blue DSS2) image of Lupus I, overlaid with ^^CO 2-1 emission from 
this work ( contours at 1, 3, 5, 7, 9Kkms~^). Annotations: da r k clou d positions (Ta- 
ble [p from [Hartley et al.l (Il986l . open triangles), IVilas-Boas et al.l (120001. fil led triangles), 
Lee fc MyersI ( 19991 . squares): embedded YSOs frorn Merin et al. ( 2008 . star-shapes); 



and 



bright stars (by name); C O 1-0 peaks from iHara et al.l ( 



1999 



diamonds). Filled squares 



denote cloud cores observed in molecular lines by iLee et al. 
^^CO map are plotted. The /, b arrows are 15' long. 



(I2OO4J ). Only objects within the 
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le'-H" 1^" 10" M" OIJ" 

Fig. 2. — Optical image of Lupus III, as in Figure [H ^^CO contours are at 2, 4, 
. . . 12Kkms~^. Bright stars are marked as crosses. 



- 39 - 




Fig. 4. — Distributions of the logarithms (base 10) of the peak intensity ratios for points 
with multiple ^^CO 2-1 (solid) and CO 4-3 (dashed) spectra, in all maps. Only ratios with 
estimated S/N>2 are included. 
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20 20 



Fig. 5. — Sample spectra toward Lupus I: each panel shows a ^^CO 2-1 spectrum (with 
position annotated). If there is a CO 4-3 spectrum toward that position (within 0.25'), it is 
shown below the ^^CO 2-1 spectrum. 
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16:13:02 
-37:53 



16:11:02 
-38:54 



O 

o 



16:08:00 
-39:10 



16:12:32 
-38:00 




16:10:05 
-38:59 






n 



20 



20 



16:11:32 
-38:08 



16:09:34 
-39:02 






O 

O 



20 



Fig. 6. — Sample spectra toward Lupus IIIN (upper row) and Lupus III (lower row). See 
Figure [5] for details. 
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o 
o 




o 



20 



Fig. 7. — Sample spectra toward Lupus IV. See Figure [5] for details. 




kig^hl Ascension U^CiQ^O) 



Fig. 8. — Maps of Lupus I in (left) ^^CO 2-1 and (right) CO 4-3, integrated from vlsr = +2 
to +10kms"^ i^CO 2-1 contours at 2, 4Kkms-^ (black) and 6, SKkms"^ (white); CO 4-3 
contours at 3, 5 (black) and 7, 9, llKkms"^ (white). The broken outline denotes the extent 
of the CO 4-3 map. 
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Fig. 9. — Maps of Lupus III and Lupus IIIN in (left) ^^CO 2-1, integrated from vlsr = +2 to 
+ 12kms^^, and (right) CO 4-3, integrated from VisR = +2 to +8kms~^ toward Lupus III 
and from +6 to +12kms^^ toward Lupus IIIN. Contours at 2, 4Kkms^^ (black) and 6, 
8. ..16Kkms-^ (white). 




Right AScSiiSieui 

Fig. 10. — Maps of Lupus IV in (left) ^^CO 2-1 and (right) CO 4-3, integrated from vlsr = 
+3 to +7kms-^ i^CO 2-1 contours at 1, 2, 3Kkms'i (black) and 4, 5, 6, TKkms"^ 
(white); CO 4-3 contours at 3Kkms~^ (black) and 5, 7Kkms~^ (white). 
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Fig. 11.— Channel maps of Lupus I in (left) ^^CO 2-1 and (right) CO 4-3, labeled with 
vlsr range in kms~^. Contours at 1 Kkms~^ (black) and 2, 3,. . . 6Kkms~^ (white). Broken 
outlines denote the extent of the CO 4-3 map. 
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Fig. 12. — Channel maps of Lupus III, as in Figure fTTl ^'^CO 2-1 contours are at 1 Kkms~^ 
(black) and 2, 4, 6, 8Kkms~^ (white); CO 4-3 contours are at 1, 2Kkms~^ (black) and 4, 
6, SKkms""*^ (white). 
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RighL Ascension (J2QO0.0) 



Fig. 13. — Channel maps of Lupus IV, as in Figure [TTl Contours at 1, 2, 3Kkms ^ (black) 
and 4, 5, OKkms"""^ (white). 
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Btght Ascension (jaijp<3.0) 



Fig. 14.— Maps of the peak of (left) ^^CO 2-1 and (right) CO 4-3 toward Lupus I. Left: 
contours at 0.5, IK (black) and 1.5, 2.5, 3.5 K (white). Right: contours at 1.5, 2K (black) 
and 2.5, 3, 3.5, 4K (white). The broken outline denotes the extent of the CO 4-3 map. 




Eight AsceiiiStdn {JSOOO Q) 

Fig. 15. — Maps of the peak toward Lupus III, as in Figure [TH contours at 1, 1.5, 2, 
2.5 K (black) and 3, 4,. . .9K (white). 
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fifiglit Ascension [jaODO 0) 



Fig. 16. — Maps of the peak toward Lupus IV, as in Figure [TH contours at 0.5, 1, 1.5 K 
(black) and 2, 2.5,. . .4.5K (white). 




Fig. 17.— Maps of the centroid velocity of (left) ^^CO 2-1 and (right) CO 4-d emission 
toward Lupus I. Contours are at 5, 5.5, 6kms~^ (black) and 6.5, 7, 7.5kms~^ (white). The 
gray scale is unsmoothed, while the contours have been smoothed with a Gaussian kernel. 
Only spectra with S/N>10 are shown. Solid outlines indicate the extent of the maps, while 
the broken outline indicates the extent of the CO 4-3 map. 
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flight Aseeiisrott (IIOWO) 



Fig. 18. — Maps of the centroid velocity of emission toward Lupus III, as in Figure [T71 
i^CO 2-1 contours are at 5, 5.5, Gkms-^ (black) and 7.5, 8, 8.5, Okms-^ (white); CO 4-3 
contours are at 4, 4.5,kms~^ (black) and 7.5, 8, 8.5, 9kms"^ (white). Only spectra with 
S/N>5 are shown. 




mgtiL A«cension {J2000.C}) 

Fig. 19. — Maps of the centroid velocity of emission toward Lupus IV, as in Figure [T71 
contours are at 4, 4.3 (black) and 4.6, 4.9, 5.2kms~"^ (white). Only spectra with S/N>10 
are shown. 



-50- 




Fig. 20.— Maps of the line widths of (left) ^''^CO 2-1 and (right) CO 4-3 toward Lupus I. 
Contours are at 2, 2.4,. . . 3.6kms^^. Gray scale is unsmoothed, while contours have been 
smoothed with a Gaussian kernel. Only spectra with S/N>10 are shown. Solid outlines 
indicate the extent of the maps, while the broken outline indicates the extent of the CO 4-3 
map. 



- 51 - 




Fig. 21. — Maps of emission line widths toward Lupus III, as in Figure [201 contours are 
at 1.5, 1.8, 2.1kms~^ (black) and 2.4, 2.7, 3kms~^ (white). Only spectra with S/N>5 are 
shown. 




Right AseensiOil CJ3000.0) 



Fig. 22. — Maps of emission line widths toward Lupus IV, as in Figure [201 contours at 1.4, 
1.6kms~^ (black) and 1.8, 2kms~^ (white). Only spectra with S/N>10 are shown. 



-52- 



LO 



CO 

o 
u 



Lu I 



o 



Q ^^^^i Vl' '■ -i I \ \ \ \ 







10 



Lu III N 



15 



10 



15 



LO 





1 




/ 


Lu III 




■ / 


/ 














/ ■■ 








/ 








/ 

■/ 








./ 

/ 






•' .' ■ • / 
• ■ ■ ■ / ■ 








■/■ 








■ /■. 








■ .■ ■■ 






^ - . ■ *. ■ . '.-y 


''.■'■/ ■ '.: '■ ' 
' ■ ' . 






■ , ;■ ■ -x'-. ■; , 
.. . ■ - O'.-' ,, J. 








" ■ V ■•■ z;*- 








L'v' V; ■'■ y^^ - ^ 









































LO 



Lu IV 



10 



15 



10 



15 



/T *di; / K-kms ^ (^^CO 2-1) 



Fig. 23. — Integrated intensities of ^^CO 2-1 and CO 4-3 for all regions. The dashed hne 
indicates a ratio of unity, while dotted hnes trace ratios (4-3/2-1) of 0.5 and 1.5. 
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Fig. 24.— Peak of ^^CO 2-1 and CO 4-3 for Lupus I and III only. The dashed hne 
indicates a ratio of unity, while dotted lines trace ratios (4-3/2-1) of 0.5, 1.5, and 3. 
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Fig. 25. — Peak of CO transitions for cores seen in C^^O 1-0. The dashed line indicates 
a line ratio of unity, while dotted lines trace ratios of 0.67, 1.5 (left), 0.33, 0.67, 1.5, 3 
(center), and 0.67, 1.5, 3 (right). Symbols denote cores in Lupus I (squares). Lup us III/III N 
(triangles), and Lupus IV (diamonds); filled symbols are C^^ O-selected cores fromlHara et al. 
( 119991 ) . and open symbols are extinction-selected cores from IVilas-Boas et al.l (120001 ). 
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Fig. 26. — Integrated intensities toward C^^O cores. The dashed hne indicates a ratio of 
unity, while dotted hnes trace ratios of 0.5 and 1.5 (left), 0.67, 1.5, 3 (center), and 3, 6, 10 
(right). Symbols are as in Figure [25l 
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Fig. 27. — Spectra toward the YSOs in Table [3|, as for Figures [5H71 



